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We present a unified description of the low temperature phase of granular metals that reveals 
a striking generality of the low temperature behaviors. Our model explains the universality of 
the low-temperature conductivity that coincides exactly with that of the homogeneously disordered 
systems and enables a straightforward derivation of low temperature characteristics of disordered 
conductors. 
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Granular metals exhibit a wealth of behaviors generic 
to strongly interacting disordered electronic systems and 
offer a unique experimental tool for studying the inter- 
play between the effects of disorder and interactions. De- 
pending on the strength of coupling between the grains 
these systems can assume either insulating- or metallic 
phases. Remarkably, metallic samples exhibit qualita- 
tively different transport properties in different temper- 
ature regimes; in particular, the low temperature phase 
appears to be similar to disordered Fermi liquids. 

The electronic transport in granular metals is governed 
by the nontrivial interplay between the diffusive intra- 
grain electron motion and grain-to-grain tunneling which 
is accompanied by sequential charging of the grains in- 
volved in the particular electron transfer process. This 
brings the notion of the Coulomb blockade, and one ex- 
pects that it is the competition of inter-grain coupling 
and electron-electron Coulomb interactions that eventu- 
ally controls transport properties of granular metals. The 
basic parameter that characterizes transport properties is 
the dimensionless tunneling conductance, g T . Depending 
on the bare tunneling conductance </i , the conductiv- 
ity can demonstrate either exponential (insulating)-, at 
< or logarithmic (metallic), at g°, tem- 

perature dependencies Q, 0, El H, experiencing metal- 
insulator transition at g§ — g°. 

The metallic phase was recently studied in Ref.Q] where 
it was shown that the low temperature dependence of the 
conductivity of granular metals coincides exactly with 
the corresponding result for the conductivity of the ho- 
mogeneously disordered samples. A question immedi- 
ately arises: is it a coincidence that two different phys- 
ical systems exhibit identical low temperature transport 
behaviors, or there is an underlying deep connection be- 
tween the two? Furthermore, do all the other physical 
quantities (specific heat, tunneling density of states, etc.) 
possess the same universality? The main result of our 
Letter is the answer to these fundamental questions. 

Generally speaking, all the universalities that one ob- 
serves in nature can usually be attributed to one or an- 
other kind of the fundamental symmetry inherent to the 



physical system in question. For example, all critical 
phenomena are described in terms of universal models 
(Ginzburg - Landau Hamiltonian) that essentially in- 
clude only the information about the large scale sym- 
metry of the order parameter corresponding to the rel- 
evant degrees of freedom. Thus in order to uncover the 
universality of inhomogeneous metals one has to seek for 
a universal description in terms of the appropriate large 
scale degrees of freedom that characterize disordered con- 
ductors. 
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FIG. 1: Schematic phase diagram for granular metals show- 
ing crossovers between three distinct phases. At small bare 
tunneling conductances, g^ < g^,, a granular metal is in 
an insulating phase at zero temperature. The metal - in- 
sulator transition in three dimensions at T = occurs at 
9t = (l/67r) ln(Ec/S), where Ec and 5 are the charging en- 
ergy and the mean energy level spacing in a single grain re- 
spectively. At large tunneling conductance, g^ > g%, the 
two different types of conductivity behavior are possible: (i) 
the high temperature phase, T > g T S is characterized by the 
logarithmic temperature dependence of the conductivity in all 
dimensions; (ii) the low temperatures phase, T < gr5 is the 
universal Granular Fermi liquid phase described in terms of 
low energy interacting diffusion modes. 
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We construct such a universal description of low tem- 
perature physical properties of granular metals building 
on the er-model introduced first for the noninteracting 
dirty metals in Refs. 0, and generalized in Ref. to 
include the interaction effects. Our approach applies 
at temperatures T < g T 6, where 5 is the mean energy 
level spacing in a single grain. The energy scale g T 5 ap- 
pears naturally as the upper energy cutoff of the effective 
model, since h/g T 8 is the mean time for the electron to 
escape from the granule ||. 

The main results of our work are as follows: Mak- 
ing use of the effective description of the granular metals 
in terms of the er-model, Eq. Ij5a(l . we show that: (i) 
All the phenomena that are described in terms of the a- 
model including interaction and localization effects, and 
all the thermodynamic quantities are universal for gran- 
ular metals, (ii) There are several physical quantities, 
which, although not directly related to the charges in the 
cr-model can, nevertheless, be found from the renormal- 
ization group equations describing the flow of charges of 
the cr-model. The important example of such a quantity 
is the tunneling density of states. We summarized our 
result in Fig. ^ where three distinct phases that one can 
identify on the basis of our approach are presented: (i) 
The universal low temperature phase, which we refer to 
as to the granular Fermi liquid, generalizes naturally the 
Fermi liquid phase of homogeneously disordered metals. 
The granular Fermi liquid phase neighbors (ii) the high 
temperature, T > g T S, metallic phase governed by the 
local single grain physics, and (iii) the insulating phase, 
where g^ 1 < g% , characterized by the activation behavior 
of the conductivity. 

Now we turn to the description of our model and the 
derivation the phase diagram in Fig. ^ We consider a 
d— dimensional array of metallic grains with the Coulomb 
interaction between electrons. The motion of electrons 
inside the grains is diffusive, and they can tunnel from 
grain to grain. We assume that in the absence of the 
Coulomb interaction, the sample would be a good metal. 
The system of weakly coupled metallic grains is described 
by the Hamiltonian 



Eq. i|la|) describes the electron-electron interaction 



H = H + H mt + ft( T i) Uj $( rj ), (la) 



where is the tunneling matrix elements (i^ = tji) 
corresponding to the points of contact and rj of z-th 
and j— th grains. The Hamiltonian Hq in Eq. (|la(l is 



^ = E/ d3r ^l [P 2 / 2 



m - + u(rj) ] fa (lb) 



with [i being the chemical potential, it describes non- 
interacting electrons scattered by random impurity po- 
tential u( r i)- The second term in the right hand side of 



Hint = -7T 2J ™i C ^ ' 



J > 



(lc) 



where CV,- is the capacitance matrix and hi = J d 3 ri fa^ fa 
is the electron number operator in the i— th grain. 

Using Eqs. Q the cr-model for granular systems can 
be derived in a usual way 0> 0] : we decouple the 
Coulomb interaction term in Eq. i|lcfl using the axillary 
fields V, average over disorder introducing Q— matrix 
field [fjQ] an d expand around the diffusive saddle point. 
The final expression for the effective low-energy action 
reads: 



S = -^£>[(e + TOi] 



* ClJ V*. 



(hi) 



2e 



2 "J" 



(2) 



Here the sums are performed over the grain indices, the 
symbol (...) means summation over the nearest neigh- 
bors, the trace is taken over spin, replica and Matsubara 
indices and e = id T . The field V in Eq. © is a vector 
in the frequency and replica spaces and the correspond- 
ing contraction is assumed: V*V = Yl u a Kj„,q Kj„,q, 
where a is the replica index and w„ = 2-nTn are the 
bosonic Matsubara frequencies. The Q— matrix in Eq. 
is the matrix in the Matsubara, spin, and replica spaces 
Q — > u , a g subject to constraint Q 2 = 1. In 

addition, each element of the Q— matrix is the quater- 
nion, i. e. it can be presented as Q — q l Ti, where 
c/i is the real vector and r, are the quaternion matri- 
ces 0. For energies e <C E T , where E T is the Thou- 
less energy, the Q— matrices in Eq. © are coordinate 
independent within each grain. In what follows we use 
Hikami-parametrization for Q— matrix |ld| 



B 

Vi - bTb 



(3) 



where the matrix B u lUn2 has non-zero elements only 
for frequencies uj ni > 0, w„ 2 < 0. Expansion of the 
Q— matrix in powers of the field B in Eq. provides a 
systematic way to take into account l/g T corrections. 

To derive the low energy a— model that provides the 
universal description of the low temperature phase, T < 
g T S, we apply the renormalization group technique to the 
effective action in Eq. J5J). We divide the field B into the 
slow and fast parts B = B s + Bf and define the fast 
part of the field B u lWn2 in a way that it exists either 
for A < u ni < A + dA or -A - dA < u n , < -A. This 
procedure (the details will be presented in the forthcom- 
ing publication) results in the following renormalization 
group equation for the tunneling conductance [2| 



dg T /dh\A = 1/(2tt(2), 



(4a) 
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where d is the dimensionality of the granular array. Inte- 
gration of Eq. (|4a|l in the energy interval (g T 8, Ec), where 
E c is the charging energy for a single grain, results in the 
following tunneling conductance 



g T =9 { r ) -(l/2Trd)ln(E c /6). 



(4b) 



For energies lower than g T S the physics is dominated by 
the distances that are much larger than the size of a 
single grain. This allows us to consider the continuum 
limit of the a— model, and arrive at the corresponding 
final expression for the action: 



4 /Tr 
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(5a) 



Here a is the grain size and the trace is taken over the 
replica indices. The renormalized diffusion coefficient D 
in Eq. (|5a|l is given by 



D = g T a 2 8 



(5b) 



with g T being the renormalized tunneling conductance 
from Eq. i|4b(l . Since the effective model l|5a|) operates 
with the Q— matrices which have only long-range degrees 
of freedom, it applies, with the appropriate charges and 
upon the high-energy renormalization, to any disordered 
metal, including a homogeneously disordered one. Thus, 
all the information about the granularity of the sample is 
hidden in the temperature independent renormalization 
of coefficients of the effective model (I5all . The conduc- 
tivity of the sample is related to the effective diffusion 
coefficient, D through the usual Einstein relation 



a = 2e 2 D(a d Sy 



(6a) 



Effective model (|5a|) together with Eq. I|4bjl for the renor- 
malized conductance naturally explains the result for the 
low temperature, T < g T 5, conductivity obtained in 
Ref. Q. The interaction correction to conductivity has 
two contributions. The first contribution is temperature 
independent and is given by 



5<Tl = -(Tq- 



1 



2ndgT 



In 



E c 



(6b) 



where cto = 2e 2 gTa 2 ~ d is the classical Drude conductivity 
of granular metals. Equation (|6b|) follows immediately 
from the renormalization of the tunneling conductance 
g T in Eq. I|4bjl and is specific to granular metals. The 
second contribution to conductivity is temperature de- 
pendent and comes from the low energy renormalization 
of the diffusion coefficient D in the effective model (15 all . 
It coincides with the corresponding correction to conduc- 
tivity obtained for homogeneously disordered metals in 



Ref. EH 
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(6c) 



where a ~ 1.83 and f3 ~ 3.13 are the numerical constants. 

Although the tunneling density of states is not directly 
related to the charges in the cr— model, it can, neverthe- 
less, be found from the renormalization group equations 
describing the flow of charges of the tr— model. Since at 
low temperatures the flow of coupling constants of the 
cr-model of granular metals is determined by the same 
renormalization group equations as in the case of homo- 
geneously disordered metals, one arrives at the important 
conclusion that the tunneling density of states has a mul- 
tiplicative structure: 



v/vq = v h v u 



(7) 



where vq is the density of states of non-interacting elec- 
trons, Vh is the contribution to the density of sates that 
comes from high energies, e > g T S, while v\ is the contri- 
bution from low energies, e < g T S which up to the proper 
renormalization of all constants coincides with the cor- 
responding result for the density of states of disordered 
homogeneous metals. As an application of Eq. J7J) let us 
consider the density of states of granular films. The low 
energy contribution, vi in Eq. J7J, is 



vi = exp 



1 1 9tS g T E% 
16g T n 2 n T n TS 3 



(8a) 



whereas the "high energy" part of the tunneling density 
of states, Vh in Eq. l|7)l. is temperature independent at 
T < g T 5 and for granular films is given by 
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(8b) 



At large tunneling conductance, g T 3> 1, the result for 
the low energy contribution, vi , in Eq. (|8a|l coincides with 
the perturbative result for the density of states of gran- 
ular metals obtained in Ref. |l|. 

Having described the effects of electron-electron inter- 
action on the transport properties of granular metals, we 
now turn to quantum- (or weak localization) corrections 
to the conductivity 12]. In the leading order in the in- 
verse tunneling conductance, l/g T , interaction and weak 
localization corrections can be considered independently. 
As usual, the weak localization correction Sa WL is defined 
by the following expression 



Sa wr , = — - 



'WL — e <?7 

7T 



C(0,q) 



d d q 
(2^p : 



(9) 



where C(0, q) is the Cooperon propagator which in 
the absence of electron-electron interaction is given by 
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with D being the effective dif- 



C(w, q) = (Dq 2 - iu 
fusion coefficient in Eq. J5bj. For 2D and ID samples 
it is important to take into account dephasing effects 
since Eq. diverges. Dephasing time may be ob- 
tained from the effective model (|5a|l straightforwardly 
using the corresponding results for homogeneously dis- 
ordered metals 13] with the proper effective diffusion 
coefficient, D = g T a 2 8. The final result for the weak 
localization corrections reads 



5a v 



c 



for granular films, and 



5a v 



2ng T 



In 



gl8_ 
T 



9t^ 
T 



1/3 



(10a) 



(10b) 



for granular wires. We notice that the quantum interfer- 
ence corrections Sa WL may be easily suppressed by ap- 
plying relatively weak magnetic field such that the main 
temperature dependence of the conductivity comes from 
electron-electron interaction effects Eqs. I|6bll and (fHcll . 

In conclusion, we have derived the low energy er-model 
that provides a universal description of the low tempera- 
ture phase of granular metals and, more generally, of any 
disordered conducting medium. This model explains a 
striking similarity of the low-temperature transport be- 
haviors of different disordered conductors as being gov- 
erned by the same long wave electronic diffusion modes. 
The proposed model enables one to derive the low en- 
ergy properties of granular metals from the corresponding 
characteristics of disordered homogeneous metals. We 
demonstrated the power of the developed approach by 



finding the density of states, the interaction and localiza- 
tion corrections to the conductivity of granular metals. 
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